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Abstract 

Using the exact path integral solution for the damped harmonic oscillator 
it is shown that in general there does not exist an exact dissipative Liouville 
operator describing the dynamics of the oscillator for arbitrary initial bath 
preparations. Exact non-stationary Liouville operators can be found only 
for particular preparations. Three physically meaningful examples are exam- 
ined. An exact new master equation is derived for thermal initial conditions. 
Second, the Liouville operator governing the time-evolution of equilibrium 
correlations is obtained. Third, factorizing initial conditions are studied. Ad- 
ditionally, one can show that there are approximate Liouville operators inde- 
pendent of the initial preparation describing the long time dynamics under 
appropriate conditions. The general form of these approximate master equa- 
tions is derived and the coefficients are determined for special cases of the 
bath spectral density including the Ohmic, Drude and weak coupling cases. 
The connection with earlier work is discussed. 
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I. INTRODUCTION 



Recently the problem of the reduced dynamics of a quantum system in contact with a 
reservoir has gained renewed interest While the foundations of quantum dissipative 

processes were layed already in the sixties 0|, this early work was mainly concerned with 
weakly damped systems and has relied on the Born and Markov approximations. In this 
conventional approach the dynamics of the dissipative quantum system is described in terms 
of quantum master or Langevin equations. Important progress beyond the limitations of the 
weak coupling approach was made in the eighties, in particular path integral techniques were 
shown to be powerful means to describe quantum dissipative systems p. This has lead to 
unexpected new results, as for instance the non-exponential decay of correlation functions 
in the low temperature range @J^, not directly available within the conventional master 
equation approach. 

The new insight gained also provided a basis for profound critique of concepts de- 

veloped in the context of the quantum master equation approach such as complete positivity 
or the quantum regression theorem. However, the results of the path integral approach were 
rarely used to derive improved master equations valid in the range of strong damping and/or 
low temperatures. In this paper we shall address this problem for the ubiquitous quantum 
dissipative system, the damped harmonic oscillator. Based on available exact results [|I^] we 
derive new generalized quantum master equations describing the exact relaxation of mean 
values and the time evolution of equilibrium correlation functions. Whenever appropriate, 
the new results will be confronted and compared with earlier findings. 

More specifically, the article is organized as follows. In the following section the micro- 
scopic model of a harmonic oscillator coupled to a harmonic bath is introduced. The exact 
time evolution of the reduced density matrix of the oscillator |10| will briefly be summarized. 
In section III we examine whether the exact time evolution of the density matrix is compat- 
ible with a generalized quantum master equation with a time dependent Liouville operator. 
It is shown that in general there is no Liouville operator independent of the initial prepa- 
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ration. In section IV we study specific types of initial preparations. Exact time-dependent 
Liouville operators describing the time evolution of thermal initial conditions and of cor- 
relation functions are derived. In the classical limit these Liouville operators are shown to 



reduce to the Adelman Fokker-Planck operator |]TT|. We also examine the case of factorizing 
initial conditions. 

In section V we investigate particular types of damping leading to a time evolution of the 
density matrix governed approximately by a time-independent Liouville operator. Earlier 
results by Talkner |^ and by Haake and Reibold are recovered. We discuss in detail the 



limit of a weakly damped oscillator and derive a generalization of the Agarwal equation [jT3| 



Additional time coarse-graining is shown to lead to the quantum-optical master equation by 
Weidlich and Haake Finally in section VI we present our conclusions. 



II. MODEL HAMILTONIAN AND EXACT TIME EVOLUTION 



The standard microscopic model |T5|-pT| for the damped harmonic oscillator starts out from 



a Hamiltonian H composed of the oscillator part 



the bath Hamiltonian 



N 2 

(2) 

and an interaction part 

N N 2 

HoK = -qT.CnXn + q'T.l^- (3) 

n=l n=l ^'"'n^n 

This model describes an oscillator with the associated classical equation of motion 

Mq{t) + t dsj{t - s)q{s) + cu^g(t) = (4) 
Jo 

where 
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7W = i|:;;^ (5) 

is the damping kernel. The dynamics of the density matrix of the entire system (oscillator 
and bath) is given by 

W{t) = exp{-tHt/n)W{0) exp{iHt/h), (6) 

from where the reduced density matrix of the oscillator is obtained by tracing out the bath 
modes 

p{t) = TinWit). (7) 

The path integral technique allows for a complete solution to this problem for a large class 
of initial states. Since the derivation is expounded in [0, we merely state the result 



Pilf^ Qf^ j dqidq[dqdq'J{qf, g}, t, g^, g-, g, g')A(gi, g-, g, q) (8) 

where A(gj, g^', g, g') is the so-called preparation function defining the initial non-equilibrium 
state Wifi) as a modification of the equilibrium state = exp(— /?if)/Tr exp(— in the 
subspace of the oscillator. This modification can be produced by operators Oj, O'^ that act 
upon the particle only leaving the reservoir coordinates unchanged 

W{Q) = Y.O,WpO'^. (9) 

We then have 

{q\Wmq) = J dqdq'\iq,q',q,q'){q\Wp\q'). (10) 

where 

A(g,g',g,g') = E(9|0,|g-)(g10;.|g'). (11) 

j 

Several relevant examples for A(g, g', g, g') are discussed in and special cases will be 



considered below, ^(g/, g^, t, g^, g^', g, g') is the propagating function describing the time 



evolution of the reduced density matrix. In terms of the transformed coordinates r = 
{q + q') /2, X = q — q' , it reads 



J {xf, Tf, t, Xi, Tj, x, r) 



1 



N{t) 



exp 



h 



j:{xf,rf,t,Xi,ri,x,r) 



where 



T^{xf, Tf, t, Xi, Ti, x,r) ^^r^ + -r^x^ 



2(g2) 2n 

+M[xfrffi{t) + Xirif2{t) - Xirff^it) - Xfuf^it)] 
+ 1-M + 2xfXiR+-{t) + x)R—{t) 



+M 



{xiC+{t) + x/Cf (t)) - ix{xiC+{t) + XfC^{t)) 



(12) 



(13) 



and 



N{t) 



27rh 



M\fs{t)\ 



2\U/2 



(14) 



The functions fi{t), R^'^{t), and Cf{t) are given in terms of the symmetrized part S{t) and 
the anti-symmetrized part A{t) — —-^G{t) of the equihbrium coordinate autocorrelation 
function C{t) = {q{t)q) = S{t) + iA(t) in the following way: 

G{t) 



flit) = h{t) 



G{ty 



City 



Mt) 

hit) = -G{t) + 

R+-{t) 
R-{t) 



G{tf 
G{t)' 



2{q')S{t)+fs{t)({qr-S'{t)) ^ 

s{t)S{t) + Mt) {{qr-s\t)) 



h [M^ {q^ 



S{t)-Mt)S{t) +{q')mt)}, 



Grit) 



(15) 
(16) 

(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



M 

T 

M 



+ f3{t)S{t) 



_M2 

sit)-ms{t) 



(23) 
(24) 



In the sequel we shall give results mainly in terms of S{t) and G{t) which both remain finite 
in the classical limit while A{t) vanishes. The Laplace transforms of G{t) and S{t) read in 
terms of the Laplace transform 7(z) of the damping kernel |p!0| , p!8 



Giz) 
S{z) 



M(3 



E 



n=-oo ^ri. ^ 



[G(^)-G(|z/„ 



Finally, 



{q') = 5(0) 



2 oo 



M(3 



n\)] ^) 



n=—oo 

M °° 



^0 + Wn\liWn\) 



are the equilibrium variances. 



III. MASTER EQUATION 



(25) 



(26) 



We now want to investigate the conditions under which the time evolution of the reduced 
density matrix can be described by a master equation 



^p{x, r, t) = C{x, r, ^, ^, t)p{x, r, t) 



(27) 



with a generally time- dependent Liouville operator C Since S is quadratic in x and r, it is 
sufficient to make the ansatz 

/ 



C = L,{t) + 



M 

T 



M \dx^ dr 
I 







( \{ 




dx 


+ ( X, r 1 




d 




) 


\ 9r J 


V 





(A 




dx 




a 


) 





X, r 



ZUt) Z„.{t) \ ( x^ 



(t) Zrr{t) 



(28) 



V) 



with complex parameter Lc{t) and complex matrices X(t), Y{t) and Z(t). Without loss 
of generality we may put Xrx{t) = Xj.r{t) and Zrx{t) = Z„.{t). Hence, there are eleven 
time-dependent functions in the ansatz (^). Using (|), ([T^) , (p!^) , and ( pSj) to write the left 
hand and right hand sides of (|27| ) in explicit form, we find that a master equation exists for 
arbitrary preparation functions provided that 15 ordinary differential equations are fulfilled. 
Three of these, obtained by comparing the coefficients of r'^^^r jTi^rf, read 

= Z„{t) + XUt)h{tf - iYrx{t)fi{t), 
Q = -iYrx{t) + 2XUt)hit). 

= xut)f^{tf- 

Inserting the solution X^^(t) = Y^.x{t) = Z„-{t) = into the 12 remaining equations, one finds 
that the conditions obtained by comparing the coefficients of 1 and VfXi imply Lc{t) = Yrr{t). 
Taking this into account, we are left with the following set of 11 equations: 



hit) = 2iXxr{t)Mtf + [Yrrit) + i;.(t)]/l(t) + 2lZxr{t) 

h{t)=2lXxr{t)U{t)h{t) 

h{t) = [Yrr{t)+2lXxr{t)h{t)]Ut) 

hit) = [Yxxit) + 2lXxr{t)fl{t)] hit) 

= -«,.,(t)/3(t)i?+-(t) + 2Xrr{t)f^{t) 

R+-{t) = [Yxx{t) + 2iXxr{t)fi{t)] R+-{t) 

+ [-2Xrrit)f,{t) - 2tXxrit)R—it)+lYxr{t)] fsit) 

R-{t) = 2 [Yxxit) + 2tXxrit)f,{t)] R-{t) 

+2ZUt) - 2iY.„{t)h{t) + 2Xrr{t)fl{t) 
Ct{t) = -2tXxr{t)Mt)C^{t) 
C^{t) = [Yxx{t) + 2zXxr{t)Mt)] C^{t) 
Ct(t) = -2iXxr{t)Ut)C:,{t) 
C^{t) = [Yxxit) + 2iXxr{t)fi{t)] C^{t) 



(29) 
(30) 
(31) 
(32) 
(33) 

(34) 

(35) 
(36) 
(37) 
(38) 
(39) 



Since only seven functions of the ansatz (^) remain to be determined, the set (p9|)-([39|) will 
be seen to have no solution in general. To demonstrate this explicitly, let us first disregard 
the equations (|36D-(p9D which stem from comparing coefficients involving the coordinates x 
and f of the preparation function. The remaining set of equations (^)-(^) has a unique 
solution and, the resulting Liouville operator can be written in the form 

To see this one first notes that (|30[) , (|3TD , and ( ^3]) give 

Xrr{t) = Yrr{t) = and X^^(t) = i/2. (41) 

The remaining four functions are then readily determined. For later convenience they are 
expressed in terms of the four functions introduced in (^Op which are given by 

G\t) - G{t)G{t) 



^q{t) = -2iZ„{t) 



G^it) - G{t)G{ty 



= .y^M = (42) 

' G^{t) - G{t)G{t) ^ ' 

Dq[t) = J^Y^rit) = lq[t){q ) ~ Jp + {q^)G{t) ^pW^W + 

n ^7 u^ (P") ^ S{t)[^p{t)X{t)+X{t)] G{t)[^,{t)Y{t) + Y{t)] 

where we have introduced 

X{t) = G{t)S{t) - G{t)S{t), 

Y{t) = G{t)S{t) - G{t)S{t). (43) 

Inserting X^rit) = i/2 and Yrr(t) = into (|36| ) -(|39| ) , it is readily seen that (|36D and ( pH]) are 
already satisfied. However, ( pT]) and (|39|) are only fulfilled provided that 

^(t)+7p(t)^(t)+7,(t)5(t) = 0, (44) 
S{t)+jp{t)S{t)+jg{t)S{t) = 0. (45) 



Differentiating (|4^) and subtracting ( ^51) we obtain 



%{t)S{t)+%{t)S{t)=0. (46) 

Now for %(t) 7^ this gives 

S{t) %{t) G{t) ^ ^ 

where the last equation follows by means of (|42|) . This implies 



S{t) = cG{t) (t>0) (48) 

where c is a real constant. Clearly, this condition, which is equivalent to Onsager's regression 
hypothesis Q], is never met exactly whatever the form of the damping kernel. To see this 
explicitly we note that the Taylor series of S{t) and G(t) start according to 

G{t) = t + 0{t^) (49) 

On the other hand, in the case 7p(t) = 0, the condition ( ^BD is only fulfilled if 7p(t) and 7g(t) 
are both independent of time. Further S{t) must be of the form 

S{t) =die-^'' + d2e-^''' (t>0), (50) 



where di and c?2 are complex constants and A1/2 = 7p/2±^y7g — 7^/4. Moreover, from (^21) 
we see that 7p(t) is constant provided that G{t) is of the form 

G{t) = cie-^i* + C2e-^2* (t > 0), (51) 

where ci and C2 are complex constants. It is now easily seen that (^) and (|51|) never hold 
exactly except for 7p = 0, which means in the absence of damping. Hence there is no exact 
master equation for the damped harmonic oscillator with a Liouville operator C independent 
of the preparation function. 

In the remainder of this work we first consider specific initial preparations for which the 
time evolution is described exactly by a time-dependent Liouville operator. We then give 
examples for approximate Liouville operators valid for particular types of damping. 
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IV. LIOUVILLE OPERATORS FOR SPECIAL INITIAL PREPARATIONS 



In the previous section we have shown that there is no exact Liouville operator which 
is independent of the preparation function. However, for certain preparations the set of 
equations (p^)-(|39|) which determine the time-dependend parameters in £ can be reduced 
allowing for an exact solution. In the sequel three types of initial states will be considered. 



A. Thermal initial condition 

Let us first consider a system that is initially in a state 

{q\Wmq')=r{q.q'){q\Wp\q') (52) 

where r(g, q') is an arbitrary function of q and q' . This initial condition allows in (^) only for 
operators Oj and O'^ that are diagonal in position space. It can be used to describe initial 
states resulting from position measurements but excludes measurements of velocities or 
variables that couple to the position and the momentum. Following Hakim and Ambegaokar 
|[T9| we call ( ^2]) a thermal initial condition. Inspection of (^) shows that the corresponding 
preparation function is given by 

A(xi, Tj, X, f) = r{xi, ri)5{x - Xi)5{f - r^). (53) 



This form of the preparation function has the consequence that only the difference of 
and (^) and likewise the difference of (|3^) and (|3^) must be fulfilled. Since (^6]) and 
are again satisfied as a consequence of (Pl), there is indeed an exact solution 

G{t)S{t) - G{t)S{t) 



G{t)S{t) - G{t)S{t) 
G{t)S{t)-G(t)S{t) 
G{t)S{t) - G{t)S{t) 



D,{t)=%{tW)-j^ (54) 
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When these coefficients are inserted into (^Oj), we obtain a time-dependent Liouville operator 
vahd for a large class of initial states. Using the relations 

1 



a; [g, ■ ], r ^ -{g, ■ }, 

d i J- X 9 ^ r 



the resulting exact master equation can be written in the form 

iM 



(55) 



Pit) 



q,l{q,pm + l{q')[p,pit)] 



h 



q.\{p.p{t)}-\{p')[q.pit)] 



(56) 



Mh 



p.\{p.pit)}-'j^{p')[q.pit)] 



Here [A, B] = AB — BA denotes the commutator and {A, B} = AB + BA the anti- 
commutator. 

To see the connection of this new master equation with earlier results, we rewrite it in 
terms of the Wigner transform of the reduced density matrix defined by 



'^{P-,(l-,t) = J dxexp (^^J^^P^ Pi^^Qjt). 



(57) 



Using the rules 



-fl-, {g,-}-2g 



hd_ 
i dq 



{p, ■ } ^ 2p 



(58) 
(59) 



we find from (E^ 



{ d 

w{p,q,t) = ig^M^iit) 



+ 



9 I 2\ 



w{p,q,t). (60) 



This is of the form of a generalized Fokker-Planck equation. For the classical harmonic 



oscillator an equation of similar form was found by Adelman based on the generalized 
classical Langevin equation 



m + ^hit) + 1^ ds-iit - s)q{s) = 



(61) 
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with a noise force ^ (t) satisfying 



m) = (62) 

and 

mm) = MkBT^{t). (63) 

The time- dependent coefficients of the Adelman equation are given by 



- Ccl(t)Ccl(t) 



^ Ccl(t)C'cl(t)-Ccl(t)Cel(t) 

C^,{t)-C,,{t)C,,{t) ^ ' 

Since in the classical limit S{t) reduces to Cci(t), and 

G{t) = -^Cciit) (65) 

the coefficients ( |5^ ) reduce to (|^). This shows that we have derived an exact quantum 
mechanical generalization of the Adelman equation. 



B. Liouville operator for the time evolution of equilibrium correlations 

In this section we investigate the time evolution of equilibrium correlation functions 

{A{t)B) = Tr (^Ae~'^'^^'BWpe'^'^^) . (66) 

where A and B are variables of the oscillator and Tr denotes the trace over the Hilbert space 
of the entire system (oscillator and bath). Inserting three partitions of unity and using 

Tr = y' cigTrR,(g| ■ (67) 

where TrR is the trace over the Hilbert space of the bath, (^) can be written as a fourfold 
integral 

{A{t)B) = J dqidq2dqsdqiA{q2, qi)B{q^, qi)P{qi, ^2, ^, ^3, Qi)- (68) 
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where A{q,q') = {q\A\q') and 

P(gi,g2,t,g3,g4) = Tir ((gi|e-^^*/^|g3)(g4|w^/3e^^*/'|g2)) • (69) 



This function contains complete information about equihbrium correlation functions. 

-P(?i5 92) ^) fjs; 94) satisfies an exact master equation. To see this we first note that an 
equilibrium correlation function may be calculated in the following way JlOl . One propagates 



the initial reduced "density matrix" Pb(0) = Bpis and takes the expectation value of A after 
time t. Hence, 

{A{t)B) = tr [ApB{t)] = J dqfdq'fA{q'f, qf)pB{qf, q'f, t) (70) 

where tr denotes the trace over the Hilbert space of the oscillator. Since the initial reduced 
"density matrix" Bpp corresponds to the "density matrix" BWp of the entire system, the 
preparation function reads 

AB(gi, q'i, g, q') = B{qi, q)S{q - g-)- (71) 

Now, using 

pBiqj, q'f, t) = J dqidq[dqdqj (g/, g}, t, g^, g-, g, q)\B{qi, g, q) (72) 

we arrive at 

{A{t)B) = J dqfdq'jdqidq[dqdq'A{qf, q'f)B{qi, q)J{qf, q'f, t, qi, q[, g, q')S{q' - q[) (73) 

which yields by comparison with (|6^) 

P{qi,q2,t,qz,qi) = j dyj{qi,q2,t,q3,y,q4,y). (74) 

Inserting the explicit form of the propagating function and after the Gaussian integration 
over y, we are left with an expression containing S{t) and G{t) only in the combination 
S{t) - ihG{t)/2M = C{t). One finds 
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P{xf,rf,t,Xi,ri) 
iM 



( iM 
'2TiN{t) ' I 



h 



M2 AC{tfN{t) 

,2 



MN{t) 



expj 

/ 



[XfVf - XiVi) + , . , . [XiVf - XfTi) 



N{t) 



N{t)C{t) 



n \C{t) 2N(t)C^(t)^ 



(75) 



where Vf = {qi + q2)/2, xj = qi - q2, Vi = {q^ + g4)/2, and Xi = q^ - q^. Further, we 
have introduced N{t) = (g^)^ — C{tY. In view of the 5 function in ([7T|) there are again 
fewer conditions that must be satisfied by the time-dependent coefficients of the Liouville 
operator. In fact, it is easily seen that 



d d 



P{Xf, Vf, t, Xi, Vi) = C{Xf, Tf, — , — , t)P{Xf, Tf, t, Xi, n), 



(76) 



with a Liouville operator C of the form (|40|) with the coefficients 

c^{t)-c{tyc{t) 



C^{t) - C{t)C{t) 
C{t)C{t) - C{t)C{t) 
C^{t) - C{t)C{t) 



(77) 



D,it)=^,it){q')- 



M2 



Hence P{xf,rf,t,Xi,ri) satisfies the exact evolution equation 

P{xf,rf,t,Xi,ri) = 



P{xf,rf,t,Xi,ri). 



(78) 



To illuminate the virtue of this equation we note that it can be used to calculate correlation 
functions in a quasi-classical manner. Introducing the double Wigner transform 

1 



P{Pf,(lf,t,Pi,qi 



dxfdXiP{xf, qf, t, Xi, q^) exp 



^{-XfPf + XiPi 



{2'KhY 

and the Wigner-Moyal transforms of A{xf,rf) and B{xi,ri) according to 



(79) 
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MPf, Qf) = J dxf exp (^-j^^fPf) M^f, Qf), (80) 
equilibrium correlations may be written as a double phase-space integral 

{A{t)B) = J dpfdqfdpidqiA{pf,qf)B{pi,qi)P{pf,qf,t,pi,qi). (81) 
This means that we can view P{pf,qf,t,pi,qi) as a quantum mechanical generalization of 



the classical joint probability. The Wigner form of the master equation ([78D, which is again 
a generalized Fokker-Planck equation, was first derived by Schramm, Jung, and Grabert 
pO| on the basis of phenomenological reasoning. Note, that the time-dependent coefficients 
are complex. However, in the classical limit they become real, since the imaginary part of 
C(t) vanishes, and the generalized Fokker-Planck equation reduces again to the Adelman 
equation. 

C. Factorizing initial preparation 



In earlier work it has been frequently assumed [Il2| , [r7i|2l| , ^ that the initial density matrix 
Wq of the entire system factorizes according to 

Wo = PoWr (82) 

where po is the density matrix of the oscillator, while Wr = ZRexp{—f3Hji) is the canonical 
density matrix of the unperturbed heat bath. Within our approach this situation cannot be 
described by a special form of the preparation function X{qi,qi,q,q'), however, factorizing 
initial states are easily gained by disregarding contributions coming from the imaginary time 
path integral (see |10| for details). The time evolution of the density matrix is then given 
by 

p{xf,rf,t) = J dridxiJFv{xf,rf,t,Xi,ri)p{xi,ri,0), (83) 



where 



Jfv(x/, r^, t, Xi, Ti) = ^[^^[^^ ^ exp[^SFv(x/, r,., t, Xi, n)] (84) 
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and 



T,FY{xf,rf,t,Xi,r,i) = M[xfrffi{t) + Xirif2{t) -Xirffs{t) -X/ri/4(t)] 



(85) 



+-M [x^iRUit) + 2xfXiR^y{t) + x}Rpy{t) . 

The index "F V" refers to Feynman and Vernon . Note that the definitions of R^^ (t) are 
now modified 



Rty{t) = h{t)[]^k,{t)-h{t)K,{t)] 



(86) 



where the fiit) are given by ([l5|)-(p!7D and 

KJt) = -J- f ds t duK'(s-u)G(s)G(u) 
M Jo Jo 

Kp{t) = j^J^ds duK\s - u)G{s)G{u) . (87) 
The function K'{t) is the real part of the real time infiuence kernel. Its Laplace transform 



is related to the Laplace transform of the damping kernel 7(z) by |10 

z 



K'iz) 



M g 



-\z^{z) 



(88) 



The integrals in (^) cannot be expressed in terms of Sit) and G(t) as it is the case for 
i?^^(t). The quadratic ansatz (^81) now leads to equations (^)-(|35|) with the functions 
R^^{i) replaced by R%^{t). There are no equations replacing (^)-([39|) since oscillator and 
bath are uncorrelated in the initial state. This set allows for a solution by the time- dependent 
parameters 

&{t)-G{t)G{t) 



G^{t) - G{t)G{t) 
G{t)G{t) - G{t)G{t) 
G^{t) -G{t)G{t) 



h 
M 

DJt) = — 



\k,{t) + ^-^k,{t) + ^,{t)K,{t) 



(89) 
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The resulting master equation is equivalent to the result by Haake and Reibold [|T2],^ who 
derived it directly from microscopic dynamics. Their formulas for Kq{t) and Kp{t) contain 
frequency integrals that may be evaluated to obtain (^) . Later, this equation was rederived 
by Hu, Paz and Zhang [^] from the path integral representation. The equivalence can most 
easily be seen using the simplified derivation given by Paz [ES . 



We mention that for factorizing initial conditions the classical limit does not yield the 
Adelman equation |^ . The generalized Fokker-Planck operator differs by terms that persist 



over times of the order of the relaxation time. This means that switching on the interaction 
with the bath at t = pathologically affects also the long time behavior of the system. Usu- 
ally, the oscillator and the bath are integral parts of the same system and the factorization 
assumption is not appropriate. 

V. LIOUVILLE OPERATORS FOR PARTICULAR TYPES OF DAMPING 

So far we have searched for exact master equations. Let us now turn to the question 
whether for particular types of damping the dynamics may be described in terms of ap- 
proximate Liouville operators valid for arbitrary preparation functions. Thus, we have to 
find circumstances under which S(t) and G(t) take the forms (|50|) and (|5TD . Using (^2]) 
we see that in this case the Liouville operator (^Oj) is time independent and the coefficients 
Dp{t) = Dp and Dq{t) = Dg read 

= (90) 



Dp = 7p 



M2 

M2' 



Note that this result is independent of the coefficients ci,C2,di,d2 in (|50D and (pT|). The 
Wigner transform of the density matrix then obeys equation ( pO| ) with time independent 
coefficients 7^ and 7^. This result is in accordance with the findings of Talkner |Q] on 
the most general form of a Liouville operator in Wigner form compatible with the correct 
equilibrium expectation values. 
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To see explicitly when conditions (0) and (|5TD hold, we first investigate the consequences 
of ([5T|). By virtue of Ehrenfest's theorem G{t) is purely classical p, since it is related to 



the response function 



xit) = 9it)^Git) (91) 



describing the mean non-equilibrium displacement in response to an applied force. Hence 



implies that the classical equation of motion is solved by a sum of two exponentials. 



This means essentially Ohmic damping. 

A. Ohmic Damping 

For strictly Ohmic damping 7(2;) = 7 and condition (51) holds exactly. We have 



G{t) = cie"^^* + C2e-^'' (92) 



where A1/2 = 7/2 rLi^uol — 7^/4 and ci/2 = =f(Ai — A2)^"'^. Thus 7p = 7 and 7^ = cug. We 
still have to examine whether S{t) fulfills (|50D. To this purpose it is useful to note that (pSf ) 
gives for the Fourier transforms 

~S{uj) = ^cot\.{^)G{u) (93) 

where 

S{uj) = S{—iu!) + S^iu), 

G{uj) =G{~tuj) -G{iuj). (94) 

The latter relations follow from the fact that S{t) is symmetric and G(t) anti-symmetric. 
Of course, (0) is just the familiar fluctuation-dissipation theorem. Now, performing the 
inverse Fourier transform we find 

S{t) = rfie-^i* + ^26"^^* - r(t), (95) 

where 



dy2 = cy,— cot (^--j, 

Tit) = ^T '^""^P^-^"^) (96) 

This means that we have to find conditions under which r{t) may be disregarded. Now 
r{t) decays at least as exp(— z/t), where u = ui = 2tx /kBT. Therefore, for temperatures 
T ^ h'y/4:7ikB, we have u ^ Re(Ai), Re(A2) and T(t) decays faster than the first two terms 
in (|95|). Hence for t ^ u^^ the last term in (|95D may be disregarded and S{t) is of the form 



However, in the strictly Ohmic case we do not have a well-defined Liouville operator 
since the sum (^6]) for (p^) is logarithmically divergent leading to associated divergences 
of the coefficients Dg and Dp in (0). To avoid this divergence we have to take the high 
frequency behavior of the damping coefficient into account which implies in realistic cases 

lim 7(z) = (97) 

z—*oo 

B. Drude Regularization 

A more realistic behavior of the damping coefficient is modeled by 

7(--) = (98) 

Z + UJd 

often referred to as Drude damping. Using (|25|) we then find 

G{t) = cie"^^* + cae-^^* + Cae'^^* (99) 



where 



Ai/2 = a±ir], As = 5 (100) 



are the solutions of 



Z -UJdZ + [uJq + ^UJd)z - UJqUJd = Q (101) 



and 
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Cl 



C2 



C3 



i a — ir] + 5 
2ri a + it] — 6^ 
i a + irj + 6 
2?7 a — i?7 — (5 ' 
2a 



(102) 



{a — ir] — 6) {a + ir] — 6) 

This result holds also in the overdamped case where r] becomes imaginary. To calculate S{t) 
we have to evaluate the inverse Fourier transform of (|9^) by contour integration. Using (|99D 
one finds 



S{t) = die-^'' + d2e~^'' + 4e 



Alt 



-As* 



m 



(103) 



where 



"j = Cj— — cot , z = 1, 2, 3 

2M I z/ / ' 



r(t) 
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E 



(104) 
(105) 



M/3„tt(A?-z.^)(Ai-z.2)(Ai-z.2)- 

From these results we see that S{t) and G{t) are of the forms (^) and (|5lD , respectively, if 
the exponentials exp(— Ast) and exp(— i/„t) decay much faster than exp(— Ai/2t). This is the 
case for 



Re(Ai),Re(A2) < 5, z/. 



(106) 



Using the Vieta relations 



2a + 5 = ojd 

2 I 2 2 / r 

+ 77^ + 2a5 = + 7cud, 



we find that Re (A 1/2) ^ 5 implies 



a -C (5. 



Further, the relations ( |107[ ) yield 



(107) 



(108) 



20 



J_ 



' 5{l + 2a/5y 



(109) 



In view of 



this gives |27 



= 2a/5 - %{a/5f + 0{{a/5f). 



(110) 



Hence, it is natural to use /ojd as small parameter to determine the roots of (|101|) . Up to 
first order in /ojd one obtains for the roots ( |100| ) 



5 = ujd ~ 2a. 



;iii) 



These relations are valid for arbitrary ratios of 7 and uoq. Also uj^jiOQ is not necessarily 
large. As a consequence of this analysis we find that only under the conditions 



7 < CUD, y 



;ii2) 



and 



;il3) 



Sii) may be approximated by the first two terms in ( |103| ). Hence, for sufficiently large 
Drude cutoff and sufficiently high temperatures, k^T ^ ^17, the oscillator dynamics can be 
described by an approximate Liouville operator with the coefficients 



2 I 2 
7g = a +7] 



7p = 2a. 



;ii4) 



This combines with (|90| ) to yield the approximate master equation 

iM 



Pit) 



h 

2i 



Mh 



ia' + v') 



q,l{q,pm+'j-{q')[p,pit)] 



q,l{p,pm-l{p')hp{t)] 



(115) 



p,^{p,pit)} 



{p')[q,Pit)] 
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first derived by Haake and Reibold |T2 



The equilibrium variances can be calculated analytically In the strictly Ohmic limit 
ujd oo, (g^) is a regular expression but (p^) diverges logarithmically. If we disregard terms 
of the order ujq/ud, i/^d, ^ I'^d-, the divergent part of (p^) is given by 

71 \ y ) 

Hence, for the classical limit it is not sufficient to have v ^ uoq, that is ksT ^ huo^, rather 
we also need '-)\\i{uj£,/v) <^ v. Thus the strictly Ohmic limit — * oo can only be taken 
after the classical limit. With this sequence of limits we obtain the classical coefficients 

The associated Liouville operator is equivalent to the classical Fokker-Planck operator of 
the Kramers equation We stress again that the results in this section remain valid for 
strong damping provided ( p.l2| ) is satisfied. 



C. Weak damping with arbitrary frequency dependence 

In this section we show that S{t) and G{t) always take the forms (^) and (|5T| ) in the 
limit of weak damping. Let us assume that the damping kernel 7(t) has a high frequency 
cutoff ujc and that its Laplace transform is an analytic function in the vicinity of —iujo. We 
introduce a typical damping strength by 

'yc= ds-f{s)cos{ujos). (117) 
Jo 

Then the weak damping condition is 

7c < cuccuc, i^. (118) 

Apart from this 7(t) is not assumed to have additional properties. 

In the limit considered one can determine the poles of G{z) from ( P^D perturbatively. To 
first order in 7c the poles are 
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Ai/2 = |±z(a;o + y) (119) 

where 

7c + «7s = 7(-^t^o)- (120) 
Performing the inverse Laplace transform of (P3| ) we find 

G{t) = (e-^^' - (121) 

with residues in zeroth order. Higher order corrections would depend on the specific form 
of 7(z) but need not be determined because the diffusion constants (|90D are independent of 
the residues. Thus, in the weak damping limit G{t) is of the form (|5l|) . From (|93D we see 
that S{lj) has the same poles as G{uj) and, in addition, poles at iun {n integer). As shown 
above, we can disregard the terms coming from the poles at ivn for times greater than the 
thermal relaxation time v^^ . Thus, S{t) is effectively of the form (0). Inserting A1/2 into 
(pop, we find to leading order in the damping strength 

jj ^ _}_ Is^O - \Vn\liW\) 

D, = ^f: 4^ = 7.^ coth (^] . (122) 

This result can be expressed in terms of the Laplace transform ( P^ of the real part of the 
influence kernel 

Dp + luo^Dq = -^K'i-zujo). (123) 

Thus, Dq and Dp are essentially given by the sine- and cosine-moments of K'{t). With 
Kc + iKg = K'{—iuJo) the master equation takes the form 

This general weak coupling master equation is given in terms of four dissipation coefficients. 
7s leads to a frequency shift and may be absorbed by renormalizing uq. 7c is the classical 
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damping coefficient. The coefficients Kc and Kg depend on the temperature. While Kc 
equals M'jc/fi' times the average energy of a quantum oscillator of frequency Uq, Kg depends 
on the specific form of ^{z) and can be calculated analytically only in certain cases. One of 
these is the Drude model. Then the moments jc/s and Kc/s are readily evaluated to read 

'c 2 I 2 

_ 7^0^D 
2 I 2 

A-„.2^|M„th(i^) (125) 
K. + -,,Re [.^,(1 + i^) - ^(1 + I 

In the strictly Ohmic limit — > oo we have 7c = 7, 7s = but 

K. = -:^\n('^) (126) 



is logarithmically divergent. The master equation (|124| ) with the coefficients ( p.25| ) can be 



compared with the well-known Agarwal equation |jT3[ 

= + ^M,^ _ « I,. _ ,„th f ^) I,, p(t)ll, (127) 



/i'2M 2 ^ /i 

which was derived with the help of projection operator techniques from the same microscopic 
model using the Born approximation in conjunction with a shot-memory approximation. As 
a main difference, we see that in Agarwal's equation the Kg term is absent. This term is 
only negligible if 

tuo < < z/. (128) 
Hence, the new master equation (|124|) is a generalization of the Agarwal equation. 



D. Connection to Lindblad theory 



The approximate time-independent Liouville operators studied above describe the dy- 
namics after the decay of fast transients. Markovian Liouville operators like these are often 
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discussed in the context of Lindblad theory |2^. This theory estabhshes the most general 



form of generators L of dissipative quantum dynamics pit) = Cp{t) preserving the positivity 
of density operators. The Lindblad master equation reads 

Pit) = -^[H,pit)] + ^j:Mt),Lt] + [h^pmt] (129) 

i 

where Li are arbitrary operators and if is a Hermitian operator. Using results by Sandulescu 
and Scutaru ||30| it is easily seen that all the above derived time-independent Liouville 



operators are not of Lindblad form. This is not too astonishing since the master equations 
derived hold only for times t > to where to is larger than an inverse cutoff frequency and u^^. 
The short time dynamics for t ^ to reduces the density matrix to a subspace where the fast 
components have decayed. The Markovian master equation holds within this subspace only, 
while Lindblad theory requires validity for any reduced density matrix. This is of course not 
necessary as has been emphasized again recently [p],|3T|. 

However, we will show that in the weak coupling limit further coarse-graining will result 
in a Lindblad operator. To this aim we first write the weak-coupling master equation ( |124|) 
in the form 

p(t) = Cp{t) = Cop{t) + iCip{t). (130) 
In terms of the usual creation and annihilation operators a\ a we have 

Co = —iido[a'^a, ■ ]. (131) 

Using the operators 

k 

where \k) are the eigenstates of a^a, Cq may be written as 

Co = ^iuJonVn- (133) 

n 

Further 
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r -^sr t 1 , 7c - ils r ^2 7c + ils ^ 2 1 
7£i = -i—[a^a, ■ ] + ^ K . " ] ^ [a , ■ ] 

+7i ([a -,01 + [a, ■ al) + 7t ([^^ " , «] + [a^ ■ «]) 



(134) 



2MujQ 

where we have introduced 

7iT 



([at . , at] + [at, ■ at]) + {[a ■ ,a] + [a,- a]) 



Kc j^Jc _Jc 

2MujQ 4 ~ 4 



coth ( ^) ± 1 



(135) 



The time evolution is formally given by 



p{t) = e"p(0). 

Now we rewrite this by a well-known operator identity 

p(t) = e^oV(O) +7 t dse'^°^'-'^Cie'^'p{<d). 
Jo 



(136) 



(137) 



For weak damping and times t <^ 7 ^ the operator e"^* in the integrand may be replaced by 
g£os_ Inserting then ( |133| ) into ( |137| ) we find 



pit) = e 



Cot 



P(0) 



(138) 



Further coarse-graining is achieved by demanding 



(139) 



Then, by performing the time integral, the off-diagonal terms [n 7^ m) are seen to be 
smaller than the diagonal terms by a factor (uot)^^. This means that in the time window 



p{t) = e^°* 1 + -ftCi p(0) 



where we have introduced the effective dissipative Liouville operator 



(140) 



:i4ii 
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The density matrix (|140|) coincides with the solution of the master equation 

p{t) = {C^ + ^C^)p{t) (142) 

for 

t < (143) 

Thus, within the time-window ujq^ <^t<^ 7^^ the two operators Cq + 'jCi and £0 + 7-^1 give 
the same dynamics. The operator Ci may be evaluated further. It is seen that only "non- 
rotating" terms containing equal numbers of creation and annihilation operators survive the 
coarse-graining in time. The resulting master equation 

p{t) = -i(a;o + y)[a^a,p(t)] +7t ([a+p(t),a] + [a^ , p{t)a]^ + 7^ ([ap(t),a+] + [a,p{t)a^]^ 

(144) 



was first derived by Weidlich and Haake |14] from a microscopic model for the damped 



motion of a single mode of the electromagnetic field in a cavity. The generator defined 
by this master equation is of Lindblad form. However, the resulting mean value equations 
violate Ehrenfest's theorem, in particular 

^ {Pit))/M. (145) 

This is due to the fact that on the coarse-grained time scale At ^ c^o^^ the variables p{t)/MujQ 
and q{t) exchange identity frequently and only a time averaged version of the mean value 
equations must be obeyed. 



VI. CONCLUSIONS 

Based on results of the path integral technique we have examined quantum master equa- 
tions for the damped harmonic oscillator. A new exact generalized master equation de- 
scribing the relaxation of initial thermal conditions was derived. This equation was shown 
to be a quantum mechanical generalization of Adelman's Fokker-Planck equation. We also 
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have given an exact Liouville operator describing the time evolution of equihbrium correla- 
tion functions which likewise reduces to the Adelman Fokker-Planck operator in the classical 
limit. The fact, that two different quantum generalizations of the Adelman operator must be 
used for the relaxation of expectation values and the regression of fluctuations is intimately 
connected with the failure of the Onsager regression hypothesis in the quantum regime. In- 
deed, the two Liouville operators (|56|) and ([76|) are only identical if S(t) is proportional to 
G{t). 

Apart from these exact results we have studied in detail the range of parameters leading to 
quantum master equations with time-independent generator. In the case of strong damping 
a time-independent Liouville operator is obtained approximately for essentially frequency- 
independent damping. However, strictly Ohmic damping is ill-behaved in the quantum case, 
and the appropriate generalization of the classical Fokker-Planck process is given by a low- 
frequency Ohmic model with high frequency cutoff such as the Drude model. This is not too 
amazing since already in the classical limit the Adelman operator becomes time-independent 
only for frequency-independent damping. 

On the other hand, in the case of weak damping the detailed frequency dependence of 
the damping coefficient is unimportant. We have derived a new generalized master equation 
valid for arbitrary weak damping as long as the sine and cosine moments of the damping 
kernel exist. The new weak-coupling master equation is more general than the well-known 
Agarwal equation. We have explained why the Liouville operator is not of Lindblad form. 
However, time coarse-graining leads to a generator of Lindblad form. The resulting coarse- 
grained master equation was found to be the Weidlich-Haake equation also known as the 
quantum optical master equation. Due to the time coarse-graining only a time-averaged 
version of the mean value equation is obeyed leading to an apparent contradiction with the 
Ehrenfest theorem. 

In summary, we have derived several new generalized master equations for the damped 
quantum oscillator for various cases of interest. In view of the new results earlier findings 
were put in proper perspective. 
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